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ABSTRACT 

Gamma-ray binaries are massive stars with compact object companions that are ob- 
served to emit most of their energy in the gamma-ray range. One of these binaries is 
known to contain a radio pulsar, PSR B1259-63. Synchrotron and inverse Compton 
emission from particles accelerated beyond the light cylinder in striped pulsar winds 
has been proposed to explain the X-ray to high energy (HE, > 100 MeV) gamma- 
ray emission from isolated pulsars. This pulsar model extends naturally to binary 
environments, where seed photons for inverse Compton scattering are provided by 
the companion star. Here, we investigate the possibility of gamma-ray emission from 
PSR B1259-63 in the framework of the striped pulsar wind model. The orbital ge- 
ometry of PSR B1259-63 is well constrained by observations and the double radio 
pulse suggests an almost orthogonal rotator so that the solid angle covered by the 
striped region is close to Air. We calculate the orbital and rotational phase-resolved 
spectral variability and light-curves to expect. We find that the recent detection by 
the Fermi/LAT of PSR B1259-63 close to periastron can be explained by a striped 
wind if the emission arises from a large range of radii IOOO77,). We constrain the 
particle density number at the light-cylinder n l ~ 7 x 10 15 m~ 3 . The re-brightening a 
month after periastron passage could be due to interaction with additional seed pho- 
tons from the trailing pulsar wind nebula. Striped winds may also be at work in the 
gamma-ray binaries LS I +61°303 and LS 5039, both of which have HE gamma-ray 
spectra reminiscent of those of pulsars and fluxes modulated on the orbital period. 
Gamma-ray pulsations are expected. Some gamma-ray binaries should be faint in HE 
gamma rays (HESS J0632+057) because the line-of-sight does not cross the striped 
wind region. 

Key words: Pulsars: general - Radiation mechanisms: non-thermal - Gamma rays: 
theory - Stars: winds, outflows 



1 INTRODUCTION 

Pulsars are thought to be highly magnetized rotating neu- 
tron stars. Many of these pulsars, previously known as ra- 
dio e mitter, have now a lso firmly been detected in gamma 
rays (|Abdo et al.ll2010al ). For some of them, the radiated 
energy is mainly released in the MeV/GeV range: there 
are called gamma-ray pulsars. Unprecedented access to this 
energy range is currently provided by the Fermi Gamma- 
ray Space Telescope Large Area Telescope. With more than 
fifty pulsars identified so far, and regular new identifications, 
compared to the handful previously known, the Fermi/LAT 
has made it possible to characterize the high energy phase- 
resolved spectra (HE, > 100 MeV) of a significant popu- 
lation of pulsars. These HE spectra are fitted by a hard 
power law, with photon index P between 1 and 2, supple- 



mented by an expone ntial cut-off at ener gies claimed to lie 
between 1 and 5 GeV (|Abdo et al.ll2010al ). Pulsed HE emis- 
sion is usually explained by some radiation processes within 
the magnetosphere of the pulsar. However, the absence of a 
super-exponential cutoff and strong gamma-ray opacities ar- 
gue against models in which the gamma-ray emission arises 
near the surfa ce of the neutron star surface ('polar cap') 
(Harding 2009), favoring models in which the HE emission 
arises from curvature radiation of particles accelerated in the 
vacuu m gaps of the outer part of the magnetosphere ('outer 
gap') (|Cheng|l2009h . 

Pulsed HE emission may also arise in the pulsar wind 
well beyond the surface of the light-cylinder (|Kirk et al.l 
2002). 3D numerical modeling of pulsar magnetospheres in 
the force- free approximation show that pulse profiles are bet- 
ter reproduced if the emission is assumed to arise in current 
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sheet s just outside of the lig ht cylinder, along the separa- 
trix (|Bai fc Spitkovskvl |2010| ). The field structure beyond 
the light-cylinder, where the rotation of the neutron star 
requires that magnetic field lines open up, tends towards a 
split monopole configuration. The misalignment of the mag- 
netic moment with respect to the rotation axis i mprints a 
spira l wave shape t o the underlying current sheet (|Coronitil 
1990; iMichel Il994l ). The resulting striped pulsar wind ex- 
pands radially at relativistic speeds with a period equal to 
that of the pulsar. Beamed emission from relativistically hot 
particles flowing in the current sheets, mostly in the first few 
stripes, leads to pulsed gamma-ray emission. In particular, 
this model has been successfully applied to reproduce ob- 
servations of the phase-resolved optical polarization of the 
Crab pulsar (iPetri fc Kir k 2005), phase-resolved spectra of 
Geminga and Vela (|Petrill2008l. ]2009h as well as light-curves 
of various Fermi pulsars (|Petrill2011T ) . 

Unlike in the outer gap models, the radiative process 
producing the HE gamma rays in the striped wind is not 
curvature radiation, which is negligible outside the light- 
cylinder, but inverse Compton (IC) scattering of ambient 
fields, typically the thermal X-ray emission from the neutron 
star itself, the cosmic microwave background or synchrotron 
self-Compton. Another possible source of seed photons, for 
pulsars in binaries, is light from the companion (let it be a 
second pulsar, a white dwarf, a low-mass or a massive star). 
The photon field orientation and density change as a func- 
tion of orbital phase, affecting IC scattering in the stripes. 
This can lead to substantial variations in the pulsed gamma- 
ray spectra along the orbit, something which is not expected 
in standard magnetospheric models where the gamma-ray 
emission arises from curvature radiation and is thus a priori 
impervious to the external photon field. 

In this paper the expected variation of pulsed gamma- 
ray emission with orbital phase is investigated using the 
specific case of PSR B1259-63, a 48 ms radio pulsar in a 
3.4 year eccentric orbit around a 10 Mq Be star (SS 2883). 
The system appears well-suited to test the model since the 
Be star has a high luminosity, ~ 10 31 W, while the or- 
bital separation d varies by more than a factor 10 between 
periastron (d ~ 0.7 a.u.) and apastron (d » 9.5 a.u.). Al- 
though pulsed emission from PSR B1259-63 has been de- 
tected in radio, only non-pulsed counterparts have been ob- 
served in X-rays and in very high energy gamma-rays (VHE, 
> 100 GeV) around periastron passage. This emission is at- 
tributed to radiation from particles accelerated at the termi- 
nation shock of the pulsar wind. Striped pulsar wind emis- 
sion would be expected to add another variable spectral 
component, peaking in HE gamma rays. Recent ly, HE emis- 
sion has been detected using the Fermi/LAT (|Abdo et al.l 
l2011al ; iTam et al.ll201ll) . A first detection coincided with 
the passage of the source around periastron (mid-December 
2011), where IC scattering from target photons of the com- 
panion star should be most efficient. The flux above 100 
MeV is ~ 10~ 7 ph cm -2 s -1 with a power-law spectrum of 
photon index ~ 2. A second detection occurred after a lull of 
about a month. The flux was about 10 times higher than in 
mid-December and the spectrum was softer (photon index 
~ 3.0). No gamma-ray pulsation were detected. 

The aim of this study is to investigate whether the HE 
gamma-ray detection of PSR B 1259-63 can result from the 
orbital phase-dependence of pulsed HE emission from the 



striped pulsar wind. The striped wind model and related 
assumptions are presented in Sj2] The specific application 
to PSR B1259-63 is described in §3] The relevance of our 
model to interpret the recent detections of PSR B 1259-63 
by Fermi/h AT and for other gamma-ray binary systems like 
LS 5039 and LS I +61°303 is discussed in Sj3]before conclud- 
ing. 



2 STRIPED PULSAR WIND MODEL 

The assumptions of the striped pulsar wi nd model, as al- 
ready discussed in several previous works l|Kirk et al.ll2002l ; 
IPetri fc Kir"k1l2005l ; [Pe"trill2009l . l201ll ) are briefly recapped in 
>i2.1l and some specificities linked to the binary environment 
are explained in 



2.1 General assumptions about the wind 

The rotation of the pulsar launches to first approximation a 
purely radially expanding wind, isotropic in space and carry- 
ing a e ± -pair density number falling off with distance to the 
center of the pulsar r faster than r~ 2 due to adiabatic cool- 
ing. These pairs are distributed according to a power law 
distribution function of index p with minimum and maxi- 
mum Lorentz factor respectively 7 m i n and 7 max such that 



dn — 72(7) dj = K e 7 p d-y 
with K e given by Eq. (8) in lPetril (2009) with 



(N - N )sech 2 (A v ij) + N 



(1) 



(2) 



This expression does not yet include adiabatic cooling be- 
cause it falls off like r~ 2 . Nevertheless, adiabatic cooling 
can be accounted for by altering the radial power law depen- 
dence of the square of the distance r in the expression for the 
particle density number. Adding an extra radial dependence, 
assuming a non accelerating wind, we find a modified power 
law such that the density number decreases as r~ 2( - p+2 ^ 3 , 
the power law exponent depends on the exact motion of the 
wind, linear acce leration would lead to another exponent 
<|Kirk et all 120021 ). Here, the maximal energy 7 max m e c of 
the pairs (in the comoving frame) sets the location of the 
cutoff in the IC gamma-ray spectrum, an important observa- 
tional constraint. We argue in §3.3.3 that the characteristic 
7 of this cutoff (or break) can be related to other parameters 
of the striped wind. No sets the background particle density 
number outside the current sheet whereas N represents the 
m axima l densi ty within the sheet. For more details, we refer 
to lPetril (|2009l ). 

The magnetic field geometry assumes a split monopole 
for which e xact analytical solutions exist, see for instance 
iBogovalovl (| 1999T ) . The current sheet separating the two re- 
gions of opposite magnetic polarity remains small compared 
to the wavelength of the striped wind, with a finite thick- 
ness parameterized by the quantity A v (defined such that its 
value is high for a thin sheet). A smooth polarity reversal 
of the toroidal component of the magnetic fie ld B v is en - 
forced in this transition layer (Eq. (5)-(6) in IPetril |2009| ) . 
Note that in the far asymptotic wind zone, the poloidal 
component of the magnetic field is negligible. We enforce 
it strictly to zero. Therefore, this asymptotic structure is a 
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simple Archimedean spiral with B r = Be 
The structure propagates radially outwards with a constant 
Lorentz factor F v and the line-of-sight is inclined by an angle 
£ with respect to the neutron star rotation axis. The sheet 
wiggles around the rotational equatorial plane between an- 
gles +x and —x, where % ls the obliquity of the pulsar, 
i.e. the angle between magnetic moment and rotation axis. 
Outside the current sheet, the plasma is cold but strongly 
magnetized whereas inside the sheet, it is hot and almost 
unmagnetized due to the zero magnetic field point (polarity 
reversal). The cold magnetized part has negligible thermal 
pressure. 

The pairs in the current sheet emit synchrotron (X- 
rays) and IC (gamma-ray) radiation. Radiative cooling of 
the particles is not taken into account but the adiabatic one 
is. However, only the first stripes contribute significantly 
since the number density of particles decreases faster than 
r~ . The emission is boosted by relativistic aberration when 
the wind expansion is directed towards the observer. Strong 
pulsed emission therefore requires that the line-of-sight cross 
the current sheet i.e. that \tt/2 — C,\ ^ X- in this special 
geometric configuration the current sheet crosses the line- 
of-sight twice per rotation of the neutron star, producing 
two pulses of boosted emission (or only one pulse in the 
less favorable case where \tt/2 — £| ~ x)- The peak-to-peak 
phase separation of the pulses is given by arccos (cot x cot £), 
which is exactly 180° for a line of sight contained in the 
equatorial plane, independent of the obliquity. The width of 
the pulses depends on the bulk Lorentz factor of the wind 
r v and on the current sheet thickness A v . The amplitude of 
the pulses is set in part by the ratio of particle density inside 
and outside the current sheet. For a detailed discussi on on 
the re lation between geometry and pulse shape, see IPetril 
<|201ll ). 

The emission extends over a range of radii starting be- 
yond the light-cylinder rz, = c P/2 tt for a pulsar of period P. 
The pulses are computed by integrating the IC emissivity 
over the first few stripes, taking into account Doppler boost- 
ing, time retardation effects and adiabatic expansion. 

The level of IC gamma-ray emission will also depend 
on the density of particles in the pulsar wind. The density 
is stronger in the striped wind region (angles between 7r/2 
and 7r/2 — x with respect to the vertical axis) because of 
the contrast N/No with the cold, strongly magnetized part 
of the wind. The maximum density for a given distance to 
the neutron star is at an angle Xi along the crest of the 
current sheet wave. The width of this maximum is tied to 
the thinness of the current sheet A v . All other things be- 
ing equal, the strongest IC emission will occur when the 
line-of-sight passes through the edge of the striped wind, in 
the region spanned around a latitude x- A single pulse in 
HE gamma rays is expected for this particular pulsar ori- 
entation. A small amount of (unpulsed) IC emission may 
also be expected when the line-of-sight does not cross the 
striped wind because of the assumed small residual density 
of particles No/r 2 away from the current sheet modulo adi- 
abatic cooling. The emission will be unp ulsed and the cold 
pairs will produce a shar p Compton line dBall fc Kirkll2000l ; 
iKhangulvan et al.1 120071 ; ICerutti et all 120081 ) or broader 
band gamma-ray emission if the pairs are assumed to 
have power law distributions (ISierpowska fc Bednarekll2005l ; 
ISierpowska-Bartosik fc Torresil2007l . 120081 ). Here, the energy 



dissipation is limited to the striped part of the wind, en- 
tailing the geometrical constraints described above for HE 
gamma-ray observability. 

2.2 Binary environment 

In a binary, the orientation of the pulsar seen by the ob- 
server does not change along the orbit (except possibly on 
super-orbital timescales due to precession of the neutron 
star). Hence, synchrotron emission from the striped wind 
stays constant. However, IC emission from up-scattered stel- 
lar photons will change with orbital phase because the angle 
between the neutron star, the companion and the observer 
changes l|Dubus et al.|[2008l ) . The calculation of the gamma- 
ray emission is changed from previous works to correctly 
follow the anisotropic scattering and to include relativistic 
aberration in the pulsar wind frame of this external source 
of seed photons. Useful analytical expressions can be derived 
for a power-law distribution of e + e~ scattering off photons 
with a blac kbody spectrum fro m a point source, both in the 
Thomson (|Dubus et all |2010|) and Klein-Nishina regimes. 
Here we use approximat e expressions for a mono ene rgetic 
photon field as given by iMoskalenko fc Strong! (|2000l ) (see 
Appendix). Relativistic boosting effects due to the pulsar 
wind motion are taken into account by Lorentz transforming 
the particle distribution function (Eq. [TJ from the co moving 
frame to the pulsar/observer frame, see IPetril (|2009l ). 



3 APPLICATION TO PSR B1259-63 

In this paragraph, we first review the essential characteris- 
tics of the PSR B1259-63 system, then show spectral fea- 
tures expected from the striped wind model before closing 
with some estimates of the relevant parameters. 



3.1 Description of the binary and pulsar 

PSR B1259-63 / SS2883 is part of the small group of gamma- 
ray binaries, systems sharing a number of observational 
properties (notably a high le vel of HE or VHE gamma-ray 
emission compared to X-rays; |Pubusll2006l ) and which have 
been conjectured to all harbor rotation-powered pulsars. 
Power-law spectra with exponential cutoffs and a flux vary- 
ing on the orbital period have been observed in two other 
gamma-ray binaries, prompting speculations as to whether 
this could be accommodated by pulsar emission (see £|-4.2p . 
PSR B1259-63 remains the only system of this group where 
a pulsar has been detected. 

PSR B1259-63 is a 47.76 ms radio pulsar in a 1236.7 day 
orbit around a 10 Mq and 6 Rq Be star (SS 288 3) with a 
temperature T, « 27, 000 K (I Johnston et al.lll994h . The dis- 
tance D to the system is about 1.5 kpc. The temperature and 
distance have been re cently revised to 27,500 -34,000 K and 
2.3 kpc respectively (|Negueruela et al]|201ll ). The orbit is 
eccentric with e = 0.87 while the inclination of the system i 
is about 36° , which implies that the distance of the pulsar to 
the Be star varies from 9.9 x 10 10 m (periastron, orbital phase 
0orb = 0) to 1.4xl0 12 m (apastron, or b = 0.5) and that the 
angle between the Be star, pulsar and observer varies be- 
tween n — i = 54° (superior conjunction, </> or b ~ 0.995) and 
n + i = 126° (inferior conjunction, or b « 0.045). Hence, the 
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stellar photon density seen by the pulsar increases by a fac- 
tor 200 at periastron, compared to apastron, and this is also 
where the effects of anisotropic IC scattering should be most 
manifest. HESS observations have detected strongly phase- 
dependent VHE emission around periastron, interpreted as 
IC emission from pairs acceler ated or randomized at the pu l- 
sar wind termination shock (jAharonian et af1l2005l , [2009). 
The first periastron passage observable by the Fermi/LAT 
occurred in mid-December 2010 and has resulted in the de- 
tection of variable HE emission (see §1). 

The measured pulse period derivative implies an energy 
loss rate for the pulsar of E f» 8 10 28 W and a magnetic 
field B ~ 3.3 10 7 T. The radio pulsations are eclipsed close 
to pe riastron by material in the circumstellar disk of the Be 
star l|Melatos et al. 1995). The pulse profile is reminiscent 
of the Crab pulse profile, with two peaks of comm ensurate 
amplitude separated by 145° l|Connors et al. I I2002T ). Pulsa- 
tions have not been detected in other wa velengths; the upper 
limit on the am plitude in X- rays is 2% (|Chernvakova et al.l 
2006). Following Ipitril l|201lf ). radio pulsations are thought 
to arise close to the polar caps, which are visible to the ob- 
server only when C ~ X- Furthermore, the double peaked 
radio pulse profile in PSR B 1259-63 suggests a nearly or- 
thogonal rotator, constraining the obliquity \ « 90°. In this 
case, the average density around the pulsar is only weakly 
dependent on latitude and pulsed HE emission is expected 
over most of the sky. The geometrical parameters of the 
striped wind model are therefore favorable and significant 
pulsed HE emission can be expected near periastron due to 
the increased seed photon density from the Be star. 

3.2 Orbital dependence of the striped wind 
emission 

From the model exposed in (JU we compute light-curves and 
spectra accounting only for photons from the companion, ne- 
glecting all other contributions. This will give a first insight 
into the orbital phase dependence of the pulse shape and 
phase-resolved spectra. When discussing the results below, 
a clear distinction must be kept between emission modu- 
lated by orbital phase due to the binarity of the system and 
modulation induced by the rotation of the neutron star and 
responsible for the pulsed emission. 

As discussed above, we expect a nearly orthogonal ro- 
tator with the line-of-sight inclination close to 90°. More 
precisely, for the following plots, we choose x = 85° and 
£ = 55°. We take a mono-energetic photon field from the 
companion with an energy of 6 eV. The wind is mildly rela- 
tivistic with T v = 10. This value for the Lorentz factor seems 
relatively small but is motivated by the fact that a similar 
Lorentz factor (r v ~ 20) was already used to satisfactorily 
fit the phase-resolve d optical polarization of the Crab pulsar 
l|Petri fc Kirk| [2005). Moreover, in our striped wind model, 
pulsed emission is assumed to be dominant very close to the 
light cylinder at the base of the relativistic outflow. This fact 
does not exclude the possibility of a bulk acceleration of the 
wind to much higher Lorentz facto r, for instance due to mag - 
netic reconnection in the stripes (|Lvubarskv fc KirkH2001f ). 
while propagating towards the termination shock. The par- 
ticle distribution function starts from 7 m i n = 10 2 up to 
/max — 10 3 with a power-law index p = 2. These values 
are justified in the discussion of <}4] 
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Figure 1. Lightcurve vs orbital phase for the pulse averaged 
high-energy emission from the striped wind above 100 MeV. The 
inset zooms in on periastron passage. 

The orbital phase light-curve is shown in Fig. [T] for 
E > 100 MeV. As expected, the flux peaks close to perias- 
tron passage (two days before). The contrast in flux is ~ 700, 
slightly more than the w 200 contrast that would be derived 
from the seed photon density because of the angle depen- 
dence of inverse Compton scattering on a point source. Note 
the asymmetrical shape with a rising time longer than the 
falling time. Significant emission occurs on a short timescale, 
with the flux greater than 10% of its peak value from 28 days 
before periastron up to 13 days after periastron. 

No variation of the pulse shape with orbital phase is 
expected because pulses are imprints of the current sheets 
that are independent of the scattering photon field, excepted 
for the peak intensity evolving with orbital phase according 
to Fig. [T] However, the shape is sensitive to the photon en- 
ergy. This variation is quantified in Fig. [2] where we plot the 
pulsed light-curves against energy from X-rays (E « 1 keV) 
to gamma-rays {E « 1 GeV, assuming electrons are accel- 
erated up these energies). At lowest energies, around a few 
keV, a non negligible off-pulse component remains visible be- 
tween both pulses whereas at highest energies, above several 
MeV, pulses sharpen with decreasing width and vanishing 
off-pulse part. We emphasize that this trend is not corre- 
lated with any orbital motion, it is an intrinsic feature of 
the pulsed radiation mechanism invoked for pulsar, entirely 
produced by relativistic beaming effects, let it be isolated or 
in binaries. In the striped wind scenario, if beaming disap- 
pears so does pulsed emission. This trend in pulse variation 
is closely related to the Doppler boosting of the spectral flux 
density from the wind frame to the observer frame. Indeed 
the Lorentz factor shapes the pulse profiles. On one side, 
if the beaming remains to weak, the DC component domi- 
nates with pulsation hardly detectable. On the other side, 
for ultrarelativistic flows, the pulse profiles become indepen- 
dent of the Lorentz factor and reflect the characteristics of 
the stripes, i.e. width and particle density number. This de- 
pendence is explained as follows. Assuming for instance a 
spectral slope of index a in intensity, the spectral flux den- 
sity F„ oc v~ a will be magnified by a factor T> 2+a where 
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log(E) (MeV) 

Figure 2. Variation of the pulsar light-curves as a function of 

energy for a given orbital phase, here we took <j> = —0.00298. A Figure 3. High-energy pulse averaged spectra of the striped wind 

smooth change with increasing contrast in pulses is recognizable for different orbital phases of the binary. The physical parameters 

when moving from lower to higher energies. Peak intensities are are estimated according to the discussion in §3.3. 

normalized to unity. 



T> — l/r v (l — /3 V ■ n) is the relativistic Doppler factor (see 
§3.3 below). The enhancement therefore strongly depends 
on this index a thus on photon energy. Note that if PSR 
B 1259-63 is an orthogonal rotator then the time lag between 
gamma-ray and radio pulses should correspond to P/4, 1/4 
of the pulsar p eriod or 90° in phase, according to Eq. (31) 
in iPetril (|201lD . 

An example of spectra for different orbital phases is 
plotted in Fig. [3] Three distinct parts corresponding to 
three different slopes are visible. First, the spectrum is hard 
(a = 0) below the low energy cut-off around a few 10 MeV 
implied by the minimum Lorentz factor of the particle power 
law distribution function. The beaming is oc V 2 . The flux 
in X-ray (around 1 keV) is two to three orders of magni- 
tude fainter than the gamma-ray flux. Most of the pulsed 
emission is radiated around 10 MeV. Second, between a few 
10 MeV and a few GeV, the classical Thomson regime ap- 
plies with a softer spectrum a = (p — l)/2, leading to more 
effective beaming oc D 2+a and sharper pulses. Finally, the 
Klein-Nishina regime corresponding to the softest part of 
the spectrum, above a few GeV, the highest values of a 
are reached, therefore a very effective beaming and strongly 
pronounced pulses. 

The characteristic frequencies of the spectrum (asso- 
ciated with the minimum electron energy and the Klein- 
Nishina transition) change slightly with orbital phase be- 
cause of the dependence of the inverse Compton cross- 
section with angle. The effect would be stronger for a 
more edge-on incli nation of the binary orbit (see Fig. 1 in 
iDubus et all (120081 )). 

3.3 Estimated striped wind luminosity 

We have reviewed the basic assumptions behind the striped 
wind model and we have shown that it leads to orbital mod- 
ulations. Here, we give some analytical estimates of the syn- 
chrotron and inverse Compton luminosity from the striped 
wind. The relevant quantities are the Lorentz factor r v , 
the magnetic field Bl and the particle density number til 
within the sheet at the light-cylinder radius tl- For PSR 
B1259-63, Bl is equal to 2.7 T at the light cylinder radius 
tl = 2.3 x 10 6 m. F v and ul will be constrained from the 
observations. 

Assuming pressure balance in the current sheet between 



the particles and the magnetic field implies 

1 / / 2 B' 2 . . 

- 7 n e m e c = - — (3) 
3 2 no 

where the 7' is the mean Lorentz factor of the pairs, n' e 
is the particle density and B' the magnetic field. Primed 
quantities refer to values measured in the frame comoving 
with the wind and the subscript L refers to values at the light 
cylinder. Introducing the magnetization parameter defined 
as 

a = = ^ = B2 (4) 

no n' e m e c 2 no Tv n e m t c 2 no Pv hl m e c 2 

we get 7' « 3<t/2. Note that although B' and n' e vary with 
distance to the pulsar, a is actually constant with radius so 
that a = ol, its value at the light cylinder radius. Therefore, 
with these assumptions, 7' is also constant with radius. 



3. 3. 1 Synchrotron flux 

The synchrotron emissivity in the comoving frame, assuming 
isotropic radiation, is 

47r 3 2 no 

assuming that the mono-energetic particles radiate only at 
their synchrotron peak (cut-off) frequency 

' — JL 12 eB ' — — ( —L\ 2 eBlj (—\ = ' (—\ 
Vs ~ 4tt 7 m e ~ 4tt V 2 ) m e F v U/""U J 

(6) 

where we have used that the magnetic field is purely toroidal. 
The synchrotron flux measured by a distant observer located 
at distance D is 

KM = 732 / 3s(v) dV^jL I V 2 j' s {v') dV (7) 

Recall that only the stripes radiate significantly. Thus the 
effective volume is much less than the total volume occupied 
by the wind, the ratio between both is symbolized by A 
which is around 0.1 and dV = Ar 2 drdQ (note that A is 
different - but related to - A^,). Integrating the above using 
the radial dependence of the synchrotron frequency (Eq. [6j 
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leads to 
F s (v) = 



Moe i S2p r 2 



9D 2 
8tt 



3 "^ e c / P A 
4 x 10" 45 Am r^ 2 Wm~ 2 Hz _1 . 



(8) 

(9) 
(10) 



where f p is a coefficient of order unity due to the integration 
over the finite solid angle covered by the striped wind fl p , 
usually less than 4n steradians, and which we collate with 
A. For a more general discussion including a power-law dis- 
tribution of part icle instead of a mono-energetic one, see 
iKirk et all (|2002l 1. 



3.3.2 Inverse Compton flux 

The inverse Compton emissivity can be approximately cal- 
culated using a similar approach and assuming the Thomson 
regime (the validity of this assumption can be verified a pos- 
teriori). If the source of seed photons (star) is a point-like 
(-R* < A) blackbody of temperature T*, the comoving ra- 
diation density of seed photons is given by 



r v (i— p v cos0*) 



0.5r 2 J m" 3 (11) 



at periastron. Here, (9* is the angle between the direction 
of the star and the direction of motion of the electron. The 
dominant emission comes from those electrons that travel 
towards the observer. Hence, the angle varies between 54° 
and 126° (see §3.1) and (1 — PvCosO*) ~ 1. In the follow- 
ing, we assume that the stellar radiation is isotropic in the 
observer frame. The striped wind emission is again assumed 
to be monochromatic at the frequency 



3.3.3 Striped wind parameters in PSR B1259-63 

The expressions above are used to constrain the parameters 
of the striped wind in PSR B1259-63. The spectrum of PSR 
B 1259-63 does not have a high enough S/N to constrain the 
presence of a n exponential cutoff during the first detection 
in December i|Abdo et al.ll2011al). when the present model 
is thought to apply (see §4). Tarn et ah! l|201l[ ) report no de- 
tection above 1.4 GeV in the first detection and a spectral 
exponential cutoff around 300 MeV during the January 2011 
flare. LS 5039 and LS I +61°303 show exponential cutoffs 
at a few GeV, in the range of the HE cutoffs (or breaks) ob- 
served in pulsars. Radiative or adiabatic losses or a complex 
distributions of pairs will readily widen the mono-energetic 
line in the simple approximation to power-law emission with 
a sharp cutoff, as observed in pulsars. We tentatively asso- 
ciate the inverse Compton frequency v c derived above with 
an emission cutoff at a few GeV, typ ical of the observed 
Fermi/LAT pulsars (|Abdo et al.ll2010al ). so that 



hv c = iVV v a L hv i , « 3 GeV. 



This implies 



10 and ul 



7 x 10 15 m" 3 . 



(17) 



(18) 



The Lorentz factor of the electrons is 7 = 3o\l2?/2 w 10 
so inverse Compton scattering on stellar photons of typi- 
cal energy e* = 2.7fcT* ~ 6 eV is indeed in the Thomson 
regime. We are effectively setting 7 max = 10 4 (in the ob- 
server frame). As an aside, the value of ul puts a limit on the 
pair multiplicity factor k in the polar caps. The Goldreich- 
Julian particle flux from both polar caps is 



27V, 



GJ 



: 47T 



£0 B ns R n 



4 x 10 32 s" 1 



(19) 



" c =3 7 "* 



-7 ,2 r v ^ = 3r v a£^ (12) 



where hv+ ~ 2.7fcT* is the peak frequency of the stellar 
blackbody in vF v . The characteristic inverse Compton fre- 
quency is independent of radius so the high energy emission 
is a line at the rest frame frequency v c = 3X>r v cr£^ w 
35 r 2 <7 2 eV with a spectral flux 

VcF c = ^ jjlf^- i<r T c 7 ' 2 < 5{u' - O AdVdv (13) 

After integration, expressing the luminosity as a function of 
T v and ul gives 

3 



VcF c 



D 2 



ulutc (T v a L ) 2 M*rf,ZemA 



1.2 x lO^r^A^Wm- 2 



(14) 
(15) 



where i cm is the range of radii over which there is emission, 
set in practice by particle cooling. In the above expression 
we have collated with A the geometrical term 



9p 



1 



47rr v 



v 3 dn 



(16) 



that appears in the integration and which is of order unity. 
Note that g p — 1 only when the striped wind covers the full 
4-7T steradians. The inverse Compton luminosity depends on 
the orbital phase through the stellar radiation density (the 
numerical value given here is for periastron). 



so the expected pair density at the light cylinder is 

2/VgJ „ 1r .10 -3 /r,r,\ 

n L w k =— w 2 x 10 Km (20 

47rr|c 

and the pair multiplicity is constrained to 

k^3x10 5 A~\ (21) 

A compara ble constraint for PSR B 125 9-63, k < 8 x 10 4 , was 
derived by IPetri fc Lvubarskvl l|2007l ) from considerations 
about magnetic dissipation at the termination shock. 

The inverse Compton luminosity at periastron directly 
depends on the adopted spectral cutoff (Eg. I14[). giving 



v c F c (3 GeV) w 1.2 x 10~ ia A ( — 



(22) 



or 4 x 10 24 W, isotropic equivalent, to compare to a pulsar 
spindown power of 8 x 10 28 W. A significant fraction of the 
power is radiated in gamma rays only if a very large range 
of radii (l eul ) are involved. At periastron, the pulsar is at a 
distance from the star equivalent to 4 x 10 ri,- Hence, a large 
fraction of the pulsar power can in principle be radiated in 
gamma rays if the striped wind propagates most of the way 
to the star. Both particle cooling in the pulsar wind and ram 
pressure from the stellar wind will act to lower the actual 
level of high energy emission. The gamma-ray luminosity 
of PSR B1259-63 at periastron represents about 5% of the 
spindown power, so ro ta IOOOt-lA -1 in order to match the 
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observed flux. Note that the emission from each pulse is 
spread over a time km/r Ec. The HE emission will remain 
pulsed only if Z C m ^ T^rL ()Kirk et al.ll2002h . which requires 
T v ^ 30A" 1/2 . 

The synchrotron spectrum is a power-law with vF v oc v. 
The synchrotron frequency is no higher than its value at the 
light cylinder (Eq. [6]) and presumably less if emission starts 
away from the light cylinder. At the light cylinder, 

v L « 2.5 x 10 11 a\ Hz, (23) 

hence, the maximum synchrotron spectral luminosity is 
(Eq.® 

v 3 F 3 (r = r L )^7x 10" 18 a|r7 2 A Wm" 2 (24) 

This must be smaller than the spindown luminosity of 
the pulsar « 3 X 1CT 12 Wm" 2 at D = 1.5 kpc, which 
yields the conservative limits ox ^ 650 T v A -1 ^ 2 and vl ^ 
10 17 rjA _1 Hz. In the above we have fixed the character- 
istic gamma-ray energy to 3 GeV, imposing T v ol ~ 10 4 . 
Therefore, the synchrotron luminosity is no larger than 
v s Fs ~ 7 x 10~ 10 r7 4 A Wm" 2 . The limit on spindown lu- 
minosity gives T v ^ 4A 1/4 and hu L < 10A~ 1/2 keV. A 
moderate Lorentz factor for the wind is enough to signifi- 
cantly reduce the synchrotron luminosity. This is consistent 
with the lack of a very hard component in the spectral en- 
ergy distribution and with the non-detection of X-ray pul- 
sations. The synchrotron flux is indeed very sensitive to the 
Lorentz factor of the wind Eq. l[8]l. Much higher value than 
those adopted in this work would imply an undetectable 
synchrotron X-ray flux. 

Taking into account the distribution in energy of pairs 
and cooling will spread the emission over a broader band. For 
instance, Thomson cooling of a mono-energetic pair injection 
results in a p = 2 distribution for the pairs and vF v oc 
v 1 / 2 . Pairs with 7 « 10 4 have an inverse Compton cooling 
timescale w 600 s at periastron so adiabatic cooling probably 
dominates over radiative cooling. Adiabatic cooling of mono- 
energetic pairs gives a p = 1 distribution and vF v oc v. The 
spectrum also has this slope below 7 m i n if the injection is 
a power-law (Fig. 4). The pulsed inverse Compton spectral 
luminosity in X-rays at 3 keV is thus expected to be from 10 3 
to 10 6 times lower than at 3 GeV, consistent with the current 
upper li mit on the X-ray pulsed fraction, 2% of 2 x 10~ 14 
W m -2 (|Chernvakova et al.ll2006l ). 



4 DISCUSSION 

The striped wind model in the environment of a massive 
star naturally leads to orbitally-modulated inverse Comp- 
ton emission in the GeV range. The model involves sev- 
eral parameters, some attributed to the pulsar wind, like 
(x, r v , N, No, A^,), supplemented with a prescribed particle 
distribution function, and some attributed to the binarity 
of the system, like star temperature, radius, orbit inclina- 
tion angle, eccentricity, periastron position and photon den- 
sity number, respectively (T*, R+, e, cj, n p h). As shown in the 
previous section, the pulsar wind parameters influence the 
observability of high energy emission, since pulsed emission 
from the striped wind can be detected only if the line-of- 
sight crosses the stripes, and the shape of this pulsed emis- 
sion (|Petrj|2oTlh . The bina ry parameters influ ence the shape 
of the orbital modulation (|Dubus et alJliooih . 



4.1 Fermi/LAT detections of PSR B1259-63 

The orbital parameters are relatively well known in PSR 
B1259-63. Pulsed emission from the striped wind can 
be detected only if the line-of-sight crosses the stripes. 
PSR B 1259-63 is close to an orthogonal rotator owing to 
its double radio-pulse structure separated by roughly half a 
period (possible only if both magnetic poles are visible). As 
a consequence, stripes should fill the whole sky and high en- 
ergy emission can be detected. Maximum inverse Compton 
flux is expected near periastron. We showed that, assuming 
the high energy emission cuts off at a few GeV like in pulsars 
or other gamma-ray binaries, the level of synchrotron and 
inverse Compton emission can be derived and is consistent 
with observations at periastron. We propose that the initial 
detection of gamma-ray emission from PSR B1259-63 is due 
to emission from the striped pulsar wind. 

Curiously, much stronger HE emission was reported 
about a month after periastron passage. The reported flux 
level « 1.5 x 10 -6 ph cm _2 s _1 is 10 times higher than dur- 
ing the first passage, c orresponding to a large fraction of th e 
spindown luminosity (|Abdo et al.l l2011al : pTam et al.ll201lh . 
Although a high gamma-ray luminosity is possible (§3.3.3), 
this is not expected to occur after periastron passage since 
both the stellar radiation density is decreasing with distance 
and the geometry for inverse Compton scattering with stel- 
lar photons becomes less favorable (increasing interaction 
angle). The increase could be unrelated to the striped wind. 
If it is, additional contributions to the seed photon field 
could be involved. 

For instance, the surrounding pulsar wind nebula 
(PWN) beyond the termination shock can contribute to the 
source of target photons. In binaries, the termination shock 
has a con e-like appearance, pointing more or less away from 
the star |Dubus et al.|[201oh . whose opening angle depends 
on the rati o of momentum fluxes from the stellar and pul- 
sar winds (Bo govalov et al.l [20081 '). The PWN emits mostly 
X-rays and gamma-rays with a luminosity bounded by the 
spindown luminosity, which is about 100 times lower than 
the luminosity of the Be star. The termination shock is typ- 
ically expected to be c l oser than the star by a factor 10 (e.g. 
iTavani fc Aronslll997l ; iDubus! 120061 ), compensating the lu- 
minosity, but the interactions are less efficient as they occur 
in the Klein-Nishina regime. However, after periastron, the 
orbital geometry is such that the pulsar gradually moves in 
between the star and the observer. Hence, the cone moves in 
between the pulsar and the observer. If the cone has a large 
opening angle, pairs moving along the line-of-sight then en- 
counter PWN photons with a progressively smaller inter- 
action angle (closer to head on) than at periastron, while 
the angle with stellar photons increases (from 54° at supe- 
rior conjunction, a week before periastron passage, to 126° 
at inferior conjunction, 55 days after periastron passage). 
As with external Compton emission in the relativistic jets 
of AGNs, the emission will be highly anisotropic. Doppler 
boosting of the PWN radiation density seen by the pairs 
moving towards the observer (in the direction of the PWN) 
and the increased upscattering rate along the line-of-sight 
can result in a much stronger gamma-ray flux. The steeper 
spectrum reported for the second detection could be due to 
the scattering occurring in the Klein-Nishina regime. 

For completeness, we also note that emission from the 
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Be disk (Ivan Soelen fc MeintiesllioTlh or thermal emission 
from the neutron star could also contribute significantly to 
the seed photon radiation density (although the scattering 
geometry is unfavorable in the latter case). We defer to fu- 
ture work an investigation of these processes and their rele- 
vance to explain the HE gamma-ray lightcurve after perias- 
tron passage. 

4.2 Evidence for a striped wind in other 
gamma-ray binaries 

Two binaries, LS 5039 and LS I +61° 303 , are established 
HE gamma-ray sources l|Abdo et al . 2009a b) and have been 
conjectured to host pulsars wit h a large spindown power 
like PSR B 1259-63 (|Dubusll2006l '). Fermi/L AT observations 
showed that the HE gamma-ray emission does not connect 
naturally to the VHE emission detected by ground based 
Cherenkov arrays above 100 GeV. In both cases, there is an 
exponential cutoff at a few GeV, typical of the HE spectra of 
a gamma-ray pulsar. The HE emission appears to be due to a 
different population of particles than in VHE regime because 
of this cutoff. A straightforward solution is to interpret the 
HE as emission coming from the vicinity of the pulsar and 
the radio/X-ray /VHE as PWN emission arising beyond the 
termination shock. 

Fermi/LAT observations showed that the HE gamma- 
ray flux from LS 5039 and LS I +61°303 is modulated on 
the orbital period, which would be difficult to explain via 
the us ual magnetospheric emission processes (|Abdo et al.l 
2009b) but is readily explained by the striped wind model 
developed in the previous sections. The HE gamma-ray 
emission should be pulsed but this is difficult to verify in 
a binary without an ephemeris for the pulsation. 

The conditions are very similar to PSR B1259-63 ex- 
cept for the smaller orbits. The radiation density at peri- 
astron is greater by a factor « 50 — 1000 and this implies 
a larger high gamma-ray luminosity. Indeed, the mean flux 
above 100 MeV is « 10" 6 ph cm" 2 s _1 . The orbital modu- 
lations are consistent with expectations from anisotropic in- 
verse Compton scattering ijAbdo et alj "2009a.b). which also 
apply to the striped wind model as discussed in the previous 
sections. 

Interestingly, the HE lightcurve of LS I +61°303 peaks 
slightly later than expected, like PSR B1259-63. As dis- 
cussed above, this might be explained if there is an addi- 
tional source of seed photons besides the star, notably the 
PWN. Changes in the source of seed photons would also 
need to be invoked to explain the long term or orbit-to- 
orbit variability in the GeV lightcurve of LS I +61°303. 
Long term variability could also be attributed to pulsar pre- 
cession. Such free precession has indeed been o bserved in 
pulsars like for PSR B1828-11 (jStairs et al.ll2000h for which 
timescales of several months have been reported. Precession 
implies changes in the orientation of the striped wind with 
respect to the orbital plane, inducing a change in the particle 
density number along the line of sight. 

A favorable orientation of the pulsar is required to in- 
voke striped wind emission, implying the number of gamma- 
ray binaries may actually be greater than currently detected. 
If the distribution of obliquities is isotropic then a ran- 
dom line-of-sight crosses the striped wind in half of the 
pulsars. The model predicts a population of sources with 



strong PWN emission (less susceptible to orientation issues) 
but weak or absent HE gamma-ray emission. Such systems 
would be most easily detected through their VHE gamma- 
ray emission by ground-based Cherenkov arrays. An exam- 
ple could be HESS J0632+057, a s ystem very strongly sus- 
pected to be a gamma-ray binary (jHinton et al . 2009) but 
that has not been detected yet by the Fermi/LAT. 

Finally, another important issue is the exact compo- 
sition of the wind. Whereas a supply of electron/positron 
is easily achieved by pair creation in the polar caps, yield- 
ing a high multiplicity factor, the only possible source for 
ions would be from the stellar crust itself. We emphasize 
that despite the much lower density of ions, they can dras- 
tically affect the dynamics of the striped wind because of 
the large separation of time and spatial scales induced by 
the large mass ratio between leptons and ions. Indeed, the 
growth rates of some instabilities can significantly deviate 
from a pure electron/positron plasma and perturb the dissi- 
pation of magnetic field lines within the stripes, modifying 
the efficiency of the bulk acceleration of the wind. However, 
there is so far no clear evidence for such an ion component. 
The possible presence of ions in the striped wind has to be 
distinguished from the hadronic interpretation of the TeV 
light-curves expected from the shocked wind. In the latter, 
interaction of the pulsar wind with the companion star disk 
material could lead to significant emission but this emission 
would be unpulsed. 



5 CONCLUSION 

We have investigated the implications of the striped pulsar 
wind model in the context of gamma-ray binaries. Inverse 
Compton upscattering of stellar photons from the compan- 
ion by pairs in the stripes generates high energy gamma-ray 
radiation that can be detected if the line-of-sight crosses the 
striped region of the wind (§2). The gamma-ray emission is 
expected to be pulsed and is also modulated along the or- 
bit because of the changing photon density and scattering 
interaction angle (§3). 

The weak HE gamma-ray emission from PSR B1259-63 
detected near periastron can be explained by striped wind 
emission. The radio pulsations suggest PSR B1259-63 is an 
orthogonal rotator so that the line-of-sight always crosses 
the stripes. The inverse Compton emission is maximum close 
to periastron. The energy of the pairs is constant with dis- 
tance from the pulsar if pressure balance between pairs and 
magnetic field is assumed in the stripes and cooling is ne- 
glected. Hence, a characteristic inverse Compton energy is 
expected that we associate with a spectral cutoff at a few 
GeV, by analogy with the spectra observed in the gamma- 
ray binaries LS 5039 and LS I +61° 303. The level of syn- 
chrotron and inverse Compton flux depends on the wind 
Lorenz factor P v , the magnetization ol and the volume oc- 
cupied by the stripes. We find that a few % of the spindown 
luminosity can be radiated in gamma-rays if the emission 
occurs over a large range of radii (§3). 

The gamma-ray emission should show a double-peaked 
pulsation at the rotation period of the neutro n star with 
the pu lse shape providing further diagnostics. lAbdo et al.l 
l|2011al ) searched for the pulsation but did not detect it (al- 
though, as they point out, PSR B1259-63 is well within the 
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range of E and Bl where pulsed gamma-ray emission is typ- 
ically detected from pulsars by the Fermi/L AT). Detecting 
the gamma-ray pulsation in this binary may be difficult to 
achieve because of the weak flux and variatio ns in timing 
solution around periastron (|Wang et al.ll2004 ). Pulsations 
may also be smeared out if the range of emitting radii is 
large and the Lorentz factor of the wind is low (§3). 

The second, much stronger HE gamma-ray detection of 
PSR B1259-63 that occurred about a month after periastron 
passage is not explained in this model without additional 
assumptions. One possibility is additional seed photons pro- 
vided by the PWN: in this case, the peak contribution arises 
after periastron when the PWN is in between the observer 
and the pulsar, a favorable configuration for inverse Comp- 
ton scattering (§4.1). We note that a similar shift in phase of 
peak HE gamma-ray emission is observed in LS I +61°303 
and could be explained similarly. 

Emission from a striped pulsar wind might solve the 
puzzle of the pulsar- like HE gamma-ray spectrum with a flux 
modulated on the orbital period in the gamma-ray binaries 
LS 5039 and LS I +61°303. The HE emission is stronger 
than in PSR B1259-63 because of the higher stellar pho- 
ton densities. Significant HE gamma-ray emission requires 
that the stripes are oriented towards us. Systems where this 
is not the case could still be detected through their radio, 
X-ray or VHE emission (attributed to the PWN) but will 
not display strong HE emission. One such system might be 
HESS J0632+057 (§4.2). 

Gamma-ray binaries may thus be providing indirect ev- 
idence for striped wind emission. If verified by the detection 
of pulsations in HE gamma rays, this would bolster the view 
that pulsar HE emission originates beyond the light cylinder. 
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APPENDIX A: ANISOTROPIC INVERSE COMPTON EMISSION FROM A POWER-LAW 
DISTRIBUTION OF ELECTRONS 

In this appendix, we give some details about the analytical expressions used to compute the spectrum arising from anisotropic 
Inverse Compton scattering by a power-law distribution of high energy leptons with index p (Eq. (0). 

The scattering rate for an anisotropic and mono-energetic target photon field with ft the direction of propagation 
of these incoming ph o tons and ei their energy, normalized to the electron rest mass energy m e c 2 , is given according to 
iMoskalenko fc Strond (|2000T ) by 



dN 
dtdE't 



(7, El, £2,C) = o 



<jt c 



£2 



£2 < 



u El (7 - £ 2 ) 2 

7 ei(l-/3-n) 
7£i (1 + /3 cosC) 
2 1 e' 1 



2-2- - + 2 +^ -^ + -+3 -^ 



l + 2e' 1 



< 7 



(Al) 

(A2) 
(A3) 

(A4) 



where 7 is the Lorentz factor of the electron, /3 its 3- velocity, £ the angle between photon direction and lepton velocity (C = 
for head-on collision) and £2 the up-scattered photon energy (also normalized to the electron rest mass energy). All useful 
quantities are measured in the observer frame except for £[ which represents the incoming photon n ormalized energy as 
measured in the electron rest frame. We checked that this formula gives identical results to Eq. (A. 5) of iDubus et al l (2008) 
in the ultra-relativistic limit 7 S> 1. In this limit, d f^ £2 is a polynomial in I/7. Integration over the particle distribution (i.e. 
over 7) can be performed analytically when the distribution is a power-law with an integer or half-integer spectral index p. 
The spectrum is given by 



dl 

dt d£2 

with 

71 = 

72 = 



(£l,£2,C) 



A" e £2 



72 dN 



dt cfe; 



■7 p d'y 



Max 



Max 



1 + W1 + 



£1 £2 (1 + cos C) 



j '"/mill 



1 + 4/1 + 



£1 £2 (1 + cos C) 



j 7 max 



(A5) 

(A6) 
(A7) 



In the case p = 2 used in our work, we find 

dN 



7 2 d'y - 



dtd£2 ' ' 12807 5 7T£l 3 (7 - £ 2 )£2 5 

(3 (60(7 ~ E2)(log( 7 ) - log( 7 - £ 2 ))7 5 + £2 (-6O7 5 + 30£ 2 7 4 + 10£ 2 V + 5£ 2 V + 3£ 2 4 7 + 2£ 2 5 )) sec 4 

107£i£ 2 (e 2 - 7) (l2(log(7 — £ 2 ) — log(7))7 4 + £ 2 (I27 3 + 6e 2 7 2 + 4£ 2 2 7 + 3e 2 3 )) sec 2 
107£i 2 £ 2 2 (e 2 - 7) (l2(log(7 - £ 2 ) - log(7))7 4 + £ 2 (I27 3 + 6e 2 7 2 - 4e 2 2 7 + 3£ 2 3 ))) 



+ 



(A8) 



